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ABSTRACT 

We present coincident observations of the Cosmic Microwave Background (CMB) from 
the Very Small Array (VSA) and Cosmic Background Imager (CBI) telescopes. The 
consistency of the full datasets is tested in the map plane and the Fourier plane, prior 
to the usual compression of CMB data into flat bandpowers. Of the three mosaics ob- 
served by each group, two are found to be in excellent agreement. In the third mosaic, 
there is a 2a discrepancy between the correlation of the data and the level expected 
from Monte Carlo simulations. This is shown to be consistent with increased phase 
calibration errors on VSA data during summer observations. We also consider the pa- 
rameter estimation method of each group. The key difference is the use of the variance 
window function in place of the bandpower window function, an approximation used 
by the VSA group. A re-evaluation of the VSA parameter estimates, using bandpower 
windows, shows that the two methods yield consistent results. 

Key words: cosmology: observations - cosmic microwave background 



1 INTRODUCTION 

Measurements of cosmic microwave background (CMB) 
anisotropics have proved invaluable in establishing the cur- 



rent ACDM model and are widely used in cosmological pa- 
rameter estimation (Spergel et al. 2003; Goldstein et al. 
2003; Rebolo et al. 2004; Readhead et al. 2004). An im- 



2 Nutan Rajguru et al. 



»»„r"*" w **** «* 

*+„ ************* **{|*f******** 
tt .fltHta***;**** 

*** ^ *+*«»^l I* 



+ ^-*** 



'* + ££*«*****' 











- "-\^\v 












~ I l : 


: 1 1 





41 *2 



Figure 1. Typical uv-coverage. Top left: VSA. Each point represents the full 1.5 GHz bandwidth. Top right: CBI. Each point represents 
a single 1 GHz channel. Bottom: Synthcsiscd beams for typical uv-coverage. 



portant check on the accuracy of CMB measurements is to 
compare the data obtained from instruments which are sub- 
ject to different systematic effects. The CBI and VSA in- 
struments have significant design differences and the com- 
parison of our data is a check that known systematic effects 
have been accurately corrected for, and that neither dataset 
is seriously contaminated by unrecognised systematic errors. 
In the CMB community, the focus has been on the compar- 
ison of power spectra. However, maps of CMB anisotropics 
also contain important information. In particular, they are 
used in tests of Gaussianity (see for example, Komatsu et al. 
2003; Aliaga et al. 2003; Savage et al. 2004), a key assump- 
tion of CMB analysis. It is therefore an important exercise 
to check the correlation between maps of CMB anisotropics. 
To this end the VSA has undertaken a programme of ob- 
serving fields previously imaged by the CBI. In this paper 
we present the results of these observations and assess their 
consistency with the CBI data. 

Key scientific results from the measurement of CMB 
anisotropics are the cosmological parameters. An essential 
ingredient in converting the flat bandpower estimates into 
cosmological parameters are the window functions which al- 
low a theoretical power spectrum to predict a flat band- 



power. Both the CBI and VSA groups use a maximum like- 
lihood method of estimating the bandpowers, with some dif- 
ferences in the implementation, but there are important dif- 
ferences in the type of window functions used for param- 
eter estimation. The CBI group computes the bandpower 
window function which fully takes into account the anti- 
correlations between neighbouring bins (Myers et al. 2003). 
The VSA group use variance windows as an approximation 
to the bandpower windows when computing parameter esti- 
mates (Rubino-Martin et al. 2003; Slosar et al. 2003; Rebolo 
et al. 2004), although anti-correlations are not accounted for. 
We assess the bias resulting from this approximation by re- 
evaluating VSA parameter estimates using bandpower win- 
dow functions. 



In section 2 we summarise the key differences between 
the VSA and CBI instruments and the implications for ob- 
serving strategies. In section 3 we describe how the differ- 
enced maps are produced and present the results of the com- 
parison. Section 4 contains our assessment of the impact of 
window functions on the VSA parameter estimation. Finally, 
in section 5 we present our conclusions. 
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Table 1. A summary of the specifications of the CBI and the VSA extended array. 





VSA 


CBI 


Observing Frequency 


33 GHz 


31.5 GHz 


Bandwidth 


1.5 GHz 


10 GHz 


Number of Channels 


1 


10 


Number of Antennas 


14 


13 


Number of Baselines 


91 


78 


Range of baseline lengths 


0.6m-2.5m 


1.0 m— 5.5 m 


t range 


R3 300 - 1500 


300 - 3500 


Primary Beam (FWHM) 


2.11° 


45.2 arcmin X 31 GHz/, 


System temperature 


fs 35 K 


» 30 K 


Mirror diameters 


0.32 m 


0.90 m 


Synthesised beam (FWHM) 


11 arcmin 


« 5 arcmin 


Flux sensitivity 


50 Jys- 1 / 2 


1.5 Jys" 1 ^ 



2 THE TELESCOPES 

The CBI is an interferometric telescope located at an al- 
titude of 5000 m in the Atacama Desert in northern Chile. 
The instrument operates in ten 1-GHz frequency bands over 
26—36 GHz. The antennas have low-noise high electron mo- 
bility transistor (HEMT) amplifiers and typical system tem- 
peratures including the CMB, ground and atmosphere are 
?s30K. The 13 Cassegrain antennas, each 0.9 m in diame- 
ter are co- mounted on a 6 m tracking platform (Padin et al. 
2002). 

The VSA is sited at the Teide Observatory, in Tener- 
ife, at an altitude of 2400 m. The VSA operates in a single 
1.5-GHz channel at a central frequency of 33 GHz. The 14 
antennas also have HEMT amplifiers and the typical system 
temperature is «35 K. In the extended array configuration, 
the VSA uses mirrors of diameter 0.322 m. The VSA horn- 
reflector antennas are mounted on a tilt table hinged along 
east- west, but each antenna individually tracks the observed 
field by rotating its horn axis perpendicularly to the table 
hinge, and wavefront coherence is maintained with an elec- 
tronic path compensator system (Watson et al. 2003). This 
is a key design difference to the CBI. 

The individual tracking of the VSA antennas allows for 
the filtering of contaminating signals. These may be celes- 
tial sources such as the Sun and Moon, or ground-spill and 
other ground based spurious signals (Watson et al. 2003). 
The fringe rates of contaminating signals differ sufficiently 
from those of the target to allow effective filtering whilst re- 
taining most of the data. For example, filtering the Sun at a 
distance of 20° from the target removes approximately 25% 
of the data. The VSA also uses a ground-shield to minimise 
ground-spill. The consequences of these design differences 
are two-fold. Firstly, the VSA is able to observe 24 hours 
a day and its extended array (as used here) can filter out 
emission from the Sun and Moon when they are as close as 
9°. The CBI is limited to observing at night and to fields 
which are more than 60° away from the Moon (Padin et al. 
2002). Secondly, the VSA is unaffected by ground-spill con- 
tamination for fields within 35° of the zenith and so is able 
to make direct images of the sky. The raw CBI data are con- 
taminated by ground-spill and this is eliminated by means of 
a differencing scheme. The method of differencing the CBI 
data used in this analysis works as follows. A lead field is 
observed followed by a trail field at the same declination but 



separated by 8 minutes in RA. The trail field visibilities are 
then subtracted from those of the lead field. This has the ef- 
fect of removing the contaminating signal which is constant 
on an 8 minute timescale, whilst preserving the statistical 
distribution of the sky Fourier modes (Padin et al. 2002). 
Both sky maps and the power spectrum are estimated from 
the differenced data. 

The CBI mounting platform allows the orientation of 
the baselines to be changed by rotating the platform about 
the optical axis. The rotation of the tracking table together 
with the broad bandwidth ensures a well-filled aperture 
plane and a circularly symmetric synthesised beam. Figure 
1 shows the typical uv-coverage and synthesised beams of 
each telescope for the observations presented in this paper, 
in the range of angular scales common to both experiments. 
The uv-coverage of the VSA is less complete and the beam 
is less circularly symmetric than that of the CBI, for these 
fields. This is partly due to the low declination of the obser- 
vations (—3.5°) which is close to the lowest elevation that 
the VSA is able to observe. The telescope is located at a 
latitude of 28° and at the higher declinations of the main 
primordial fields, uv-coverage is significantly better (Taylor 
et al. 2003; Dickinson et al. 2004). Unlike the CBI, the VSA 
mounting table does not allow for rotation about the optical 
axis. Together with the lower bandwidth this limits the cov- 
erage of the aperture plane, although minimises the number 
of redundant baselines at a given frequency. A summary of 
the specification of the two telescopes is shown in Table 1. 



3 OBSERVATIONS 

The CBI data used in this comparison are the mosaics 02H, 
14H and 20H and first reported in Pearson et al. (2003). The 
CBI data also include two deep field observations which fall 
within the 14H and 20H mosaics and are described in Mason 
et al. (2003). Both CBI mosaiced and deep field data were 
collected during the period January - December 2000. Each 
mosaic consists of 42 differenced fields. 

The VSA observations were carried out at a later epoch. 
The data were made between May 2002 and May 2004 with 
the telescope in the extended array configuration. The larger 
primary beam size of the extended array mirrors allows each 
mosaic to be covered in three pointings. Table 2 shows the 
coordinates and effective integration times for the VSA ob- 
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Table 2. Celestial coordinates for the VSA observations of the 
CBI 02H, 14H and 20H mosaics. The effective integration time is 
calculated after flagging and filtering of the data. 



Field RA (J2000) DEC (J2000) t int (hrs) 



VSA-02H-a 


02 


44 


24.00 


-03 


30 


00.0 


169.9 


VSA-02H-b 


02 


50 


00.00 


-03 


30 


00.0 


86.2 


VSA-02H-C 


02 


55 


36.00 


-03 


30 


00.0 


49.7 


VSA-14H-a 


14 


44 


24.00 


-03 


30 


00.0 


37.1 


VSA-14H-b 


14 


50 


00.00 


-03 


30 


00.0 


73.3 


VSA-14H-C 


14 


55 


36.00 


-03 


30 


00.0 


44.9 


VSA-20H-a 


20 


44 


24.00 


-03 


30 


00.0 


154.2 


VSA-20H-b 


20 


50 


00.00 


-03 


30 


00.0 


94.2 


VSA-20H-C 


20 


55 


36.00 


-03 


30 


00.0 


74.4 



servations. The VSA data reduction and calibration proce- 
dure are described in Dickinson et al. (2004) and references 
therein. 

3.1 Foreground contamination 

At frequencies of 26-36 GHz, the dominant cosmological con- 
tamination to CMB observations comes from galactic fore- 
grounds and extragalactic radio sources. The diffuse galactic 
foregrounds include both bremsstrahlung and synchrotron 
emission. However, this emission is concentrated in the 
galactic plane and contamination may be avoided by observ- 
ing at high galactic latitudes. The observations presented 
here are at galactic latitudes > 20°. In addition, the data 
are insensitive to the large angular scales where galactic con- 
tamination is significant. Therefore, the main contaminant 
is likely to be extragalactic point sources. 

The standard VSA source-subtraction strategy involves 
an initial survey with the Ryle Telescope (RT) at 15 GHz 
(Waldram et al. 2003). Sources identified by the RT are then 
monitored with the source subtractor at 33 GHz. The obser- 
vations are carried out simultaneously with the CMB field 
observations to take account of the variability of the sources. 
A statistical correction is also applied to the power spectrum 
to remove the small effect of the remaining, fainter sources. 
As the RT is located in Cambridge at a latitude of +52°, we 
were unable to survey fields at the low declinations of the 
CBI observations. For this reason, a more limited level of 
source subtraction was implemented. This involved the sub- 
traction from the data of both groups, the fluxes obtained 
by the CBI group from observations at 31 GHz with the 
40-metre telescope at the Owens Valley Radio Observatory 
(OVRO). These observations were carried out, simultane- 
ously, in so far as was possible, with the CBI observations. 
As noted above, the VSA observations were carried out a 
later epoch and no further source observations were carried 
out to account for the variability of sources. 

The typical sensitivity achieved in the OVRO data was 
2mJy (rms) which allows for a completeness estimate of 
90% at S 3 i > 16mJy (Mason et al. 2003). To achieve high 
i measurements of the power spectrum, a deeper level of 
discrete source subtraction is required. To achieve this the 
usual approach of the CBI group is to employ the strat- 
egy of constraint matrices to 'project out' sources at known 
positions but with unknown fluxes (Bond et al. 1998). All 



sources >3.4mJy in the 1.4 GHz NVSS catalogue are pro- 
jected out of the CBI data. A statistical correction based on 
the CBI source count is then applied for sources < 3.4 mjy 
at 1.4 GHz (Mason et al. 2003). For this investigation, the 
OVRO fluxes have been subtracted from the CBI data but 
the constraint matrix strategy has not been employed. This 
is due to the limited resolution of the data included in the 
comparison and to avoid removing a significant fraction of 
the data, which is unnecessary for this analysis. 

The different methods usually employed to remove the 
effects of contaminating sources are a key distinction in the 
analysis of the two groups. Although we are unable to make 
a direct comparison of these strategies, due to the positions 
of the fields, the 33 GHz source counts estimated from the 
main VSA source monitoring programme have been used to 
evaluate the CBI source subtraction strategy. Based on VSA 
source counts Cleary et al. (2004) find the residual correction 
due to sources below the detection threshold of 3.4 mjy at 
1.4 GHz to be 0.03 Jy 2 sr _1 which is consistent with the CBI 
group estimate of 0.08 ± 0.04Jy 2 sr~\ 

3.2 Maps 

As a consequence of the need to difference the CBI visibil- 
ities, maps from the two telescopes may only be compared 
if the differencing scheme is also implemented on the VSA 
data. Further to the usual data reduction, the VSA data 
were processed is several ways to enable this. 

The VSA has a much larger primary beam than the 
CBI and covers the same area of sky in fewer pointings. To 
produce data equivalent to each CBI lead and trail field, the 
following procedure was implemented. The first step was to 
shift the VSA field centres to each of the CBI field centres. 
The complex visibility measured by an interferometer is de- 
fined as: 

V(u) = J d 2 x A(x) I(x) e - 2 ™ u - x + N(u) (1) 

where A(x) is the primary beam, I(x) is sky intensity, 
u = (u, v) is the baseline length, measured in units of the 
wavelength and N(u) is the instrumental noise (Thompson 
et al. 2001). The field shift was achieved by rotating the 
phase of the VSA visibilities so that the direction cosines, 
x — (Ax, Ay), were defined with respect to a new phase 
centre: 

Ay = sin S cos<5o — cos<5 sin<5o cos(a — ao) 
Ax = cos 5 sin(a — ao) 

where a, S and ao, So are the right ascension and declination 
of the VSA and CBI field centres respectively. 

At this stage, before differencing can be carried out, 
the larger primary beam of the VSA must be taken into ac- 
count. The effect of observing a limited area of sky is to con- 
volve the sky Fourier modes with the aperture function (the 
Fourier transform of the primary beam) . To mimic observa- 
tions with the CBI beam which has a FWHM of 45.2 arcmin- 
utes x (31 GHz/v), the VSA visibilities were convolved with 
the CBI aperture function. This was modelled as a Gaussian 
with a cut-off at 0.45 m (the outer radius of the antenna) and 
a central region with no illumination at r < 0.0774 m (corre- 
sponding to blockage by the secondary mirror). The 84 lead 
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Figure 2. The central region of the 42-field mosaics of differenced maps. Each map covers an area of 2.4° x 2.4°. The RA scale refers 
to the position of the lead field. Left: VSA data Right: CBI data. Top: 02H mosaic Centre: 14H mosaic Bottom: 20H mosaic 
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Figure 3. Sensitivity maps for each mosaic shown in 2. Each map 
ISmJybeam" 1 . The RA scale refers to the position of the lead field, 
mosaic Bottom: 20H mosaic 



covers an area of 2.4° x 2.4° and is blanked at a threshold of 
Left: VSA data Right: CBI data. Top: 02H mosaic Centre: 14H 
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and trail fields for each mosaic were produced by convolv- 
ing the phase rotated data with the CBI aperture function. 
Corresponding lead and trail fields were then differenced. 

It is important to note that matching the uv-coverage 
of the two interferometers is critical in making maps of the 
CMB. This is a consequence of sampling a random Gaus- 
sian field. Not only must the same range of angular scales 
be used, but they must sample the same region of uv-space. 
Figure 1 illustrates the mismatch in uv-coverage between 
the VSA and CBI in the range of common uv-scales. To 
overcome this, the data were binned and re-weighted. As 
the data were obtained from a convolution of the sky fourier 
modes with the CBI aperture function, which has a FWHM 
of 67A, a cellsize of 17A was chosen to ensure that the aper- 
ture function was more than adequately sampled according 
to the Nyquist sampling theorem. Data cells with a match 
in both sets were assigned a weight of the geometric mean of 
their individual weights. Data cells adjacent to, but without 
a direct match, were downweighted by a factor of two. Cells 
without a direct or adjacent match were assigned a weight 
of zero. 

Maps were then produced from the re-weighted visibili- 
ties using the AIPS package and are shown in Figure 2. The 
CBI data have been standardised to 33 GHz with the as- 
sumption of a blackbody spectrum. Both datasets have been 
corrected for the primary beam. In the case of the VSA data, 
the data were corrected for an effective primary beam where 
the effective beamsize, <7 c ff, is given by: 

1 11 

°eff a C °V 

where ac and <rv are the CBI and VSA primary beam sizes 
respectively. The effective beam size is smaller than the CBI 
beam size by 5%. This is because the sky signal measured 
by the VSA data has been multiplied by the CBI and VSA 
primary beams. 

Figure 3 shows the sensitivity maps for each mosaic. 
The minimum noise in the VSA 02H, 14H and 20H mosaics 
is 7.29 mjy beam -1 , 7.44 mjy beam -1 and 6.03 mjy beam -1 
respectively. The corresponding noise levels for the 
CBI mosaics are 3.88 mjy beam -1 , 1.82 mjy beam -1 and 
1.08 mjy beam -1 . The noise levels in the CBI maps are 
approximately a factor of 2 lower than those of the VSA 
maps, with a wider gap in the region of the CBI deep fields 
which are clearly visible in the sensitivity maps. The blank- 
ing threshold of 15mJybeam -1 is for illustrative purposes 
only. In calculating the map correlations, a cut was applied 
when the power in the primary beam reached 1/e of the 
maximum level. 

There are a number of statistics which may be used 
to quantify the consistency between datasets. The measure 
used here for both the map plane and uv-plane tests is the 
product-moment correlation coefficient (e.g. Barlow 1989). 

_xy — xy 
<j x a y 

In both planes a weighting scheme was used. In cor- 
relating the maps, each pixel was weighted by the power of 
the primary beam at the centre of the pixel. In the uv-plane, 
each cell was weighted by the inverse of the noise squared. As 



Table 3. Correlation of the 02H and 14H mosaics. The expected 
correlation is in excellent agreement with the actual correlation 
of the 02H and 14H mosaics. 





Expected 


02H 


Map plane 


0.23 ±0.02 


0.25 ±0.03 


uv-planc 


0.12 ±0.02 


0.14 ±0.02 




Expected 


14H 


Map plane 


0.23 ±0.02 


0.24 ±0.03 


uv-plane 


0.12 ±0.02 


0.13 ±0.02 



Table 4. Correlation of the 20H mosaic. The expected correlation 
in the 20H mosaic given the greater weight of data in these fields 
is shown in the first column, and in the second column this is 
modified by the inclusion of the VSA summer phase errors. The 
actual correlation, in column three, is consistent with the level 
expected given the phase calibration errors. 



Expected Expected 20H 
(inch phase errors) 

Map plane 0.27 ± 0.02 0.23 ± 0.02 0.22 ± 0.03 

uv-plane 0.14 ±0.02 0.11 ±0.02 0.10 ±0.02 



the sky is real, the underlying fourier modes a(u) = o*(— u). 
Consequently, correlations exist between visibilities that lie 
on opposite sides of the uv-plane. To account for this, the 
conjugate symmetry of the visibilities was used to reflect the 
data into one half of the plane and the correlation was car- 
ried out only in this region. The real and imaginary parts of 
the visibilities were treated independently. 

One of the disadvantages usually cited with this mea- 
sure of correlation is the difficulty of interpretation. This 
has been overcome by using Monte-Carlo simulations of the 
data to predict the expected correlation. The input model 
for these simulations was a standard ACDM model with 
noise levels and uv-coverage appropriate to the actual ob- 
servations. The method of producing the VSA differenced 
data from the simulations was carried out in the same man- 
ner as for the actual observations. 

The correlation of the 02H and 14H mosaics is shown in 
Table 3. There is excellent agreement between the observed 
and expected correlations in the 02H and 14H mosaics in 
both the map plane and the uv-plane. 

The expected correlation of the 20H mosaics is higher 
than for the 02H and 14H mosaics, which reflects the greater 
signal-to-noise in these observations. The actual correlation 
however, lies 2<r away from the expected level - see Table 
4. In contrast to the first two mosaics, the VSA data in the 
20H mosaic were largely collected during the summer period 
and between the hours of 8am and 6pm. In this period, the 
accuracy of the phase calibration is known to deteriorate 
and typical phase errors are some 20°. The source of these 
errors is thought to be due to a slight warping of the tilt 
table in the heat. All of the VSA data in the 02H and 14H 
mosaics were made outside of the summer daytime period, 
when typical phase errors are around 3-4°. 

The VSA phase errors are estimated from the change 
in phase calibration factors over the course of a day. Sev- 
eral calibration observations are carried out each day but 
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Figure 4. The power spectrum of the 02H, 14H and 20H mo- 
saics from the undifferenced VSA data (open circles) and the dif- 
ferenced CBI data (solid circles). The concordance model which 
best fits WMAP data is shown for comparison. 



^ 0.01 - 
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Table 5. The CMB band-powers (T{*£(£ + l)C e /2ir [lJ.K 2 ]) for the 
joint mosaic power spectra shown in figure 4. £ e g is the centroid 
of the bandpower window function. 



Bin 


f-range 


VSA 


/VSA 
efT 


CBI 


err 


1 


- 


-400 


3730 ± 669 


306 


3455 ± 973 


311 


2 


400 


- 586 


1798 ± 372 


487 


1993 ± 432 


490 


3 


586 


- 772 


1858 ± 386 


679 


1955 ± 367 


673 


1 


772 


- 957 


2363 ± 490 


861 


2396 ± 405 


858 


5 


957 


- 1143 


1382 ± 580 


1044 


1099 ± 244 


1045 


6 


1143 


- 1329 


1709 ± 907 


1232 


1873 ± 309 


1226 


7 


1329 


- 1515 


15 ± 817 


1416 


1055 ± 251 


1434 



there may be significant change in the phase calibration be- 
tween observations. The calibrators are either point sources 
or resolved sources for which a model is known. The phase 
errors have a weak dependence on baseline length but show 
a high degree of repeatability. The expected correlation is, of 
course, revised downwards by modelling the summer phase 
errors in the simulations. The phase errors fully account for 
the reduced correlation seen in the 20H mosaics. 

In principle, the power spectrum is insensitive to the 
phases of the visibilities, although this depends upon the 
level of correlation of the phase errors. The phase errors on 
the VSA data are found to be highly correlated. Although 
typical summer phase errors are some 20° , the rms about the 
mean phase error for each baseline is only 2-3°. A further 
consideration is the number of baselines which contribute to 
each uv-cell. If a large number of cells have contributions 
from several baselines then the phase errors will be less cor- 
related and it is possible that this could affect the power 
spectrum estimate. However, simulations of VSA shallow 
field observations which include the effect of the crossing 
uv-tracks and summer phase errors show that this has a 
negligible effect on the power spectrum. However, as a fur- 
ther safeguard the VSA group discard the worst affected 
portion of summer data. 

Figure 4 shows the joint power spectrum for the 02H, 
14H and 20H mosaics in the range of scales common to both 
instruments. The VSA power spectrum was estimated using 
the direct VSA observations, prior to the implementation of 
the differencing and reweighting scheme. There is excellent 
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Figure 5. top: VSA variance window functions, bottom: VSA 
bandpower window functions. Alternate curves are solid and dot- 
ted lines for clarity. 



agreement in the first six bins. As the contribution of point 
sources increases as £ 2 , the discrepancy in the final bin may 
be due to the variability of sources, as the data were col- 
lected at different epochs. Table 5 shows the binning and 
bandpowers for the power spectra. 



4 WINDOW FUNCTIONS 

It is standard practice for experiments observing a limited 
sky area to assume a flat bandpower when estimating the 
CMB power spectrum (Kuo et al. 2002; Hobson & Maisinger 
2002; Myers et al. 2003). However, the theoretical power 
spectra are not flat and in estimating the cosmological pa- 
rameters it is necessary to define a window function which al- 
lows a theoretical prediction of the flat bandpower qs which 
may be compared to the experimental values. The theoreti- 
cal prediction is the expectation value of qs for a given input 
model: 
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Ce (a p ) 



(2) 
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To estimate the bandpowers, the CBI group use the second 
order Taylor expansion of the log likelihood function around 
the maximum likelihood bandpowers. A quadratic approx- 
imation for the change in the bandpower, Sqs, is used it- 
eratively, to move towards the maximum, with the second 
derivative of the log likelihood function replaced by its ex- 
pectation value, the Fisher matrix. The expectation value of 
the bandpowers is then given by 



(3) 



where C is the covariance matrix with independent con- 
tributions from all signal sources: the CMB signal C , the 
instrumental noise C N and the foreground signals C src and 
C rcs . C|, is the CMB signal from each band (Myers et al. 
2003). Since 



B 



(4) 



the bandpower window functions can be calculated, once the 
maximum likelihood bandpowers have been obtained, from 



(c-^c- 1 ) 



dC b 



dC e 



(5) 



where F is the Fisher matrix of the fine-binned Ci estimates 
(Knox 1999). 

The VSA group also use a maximum likelihood method 
to estimate the bandpowers but do not use the approxi- 
mations described above. Instead the exact likelihood func- 
tion is calculated as a function of each bandpower through 
the maximum-likelihood point, with some speed-up achieved 
by the use of the signal-to-noise eigenbasis (Hobson & 
Maisinger 2002). For parameter estimates the VSA group 
do not compute the bandpower window functions but use 
as an approximation the variance window function which is 
rapidly constructed from the overlap integrals of the aper- 
ture functions for pairs of visibilities and is defined as 



l Jo 



|S(£,0)| 2 d0, 



where 



27r|u| and 



S(u) = ^ w fc i eff (u - v fc ) 



(6) 



(7) 



where u>k is the noise weighting of the visibilities and A e ^ 
is the effective aperture function (Scott et al. 2003) 

Knox (1999) has noted that the bandpower window 
function is distinct from the variance window function where 
the signal or noise is correlated. Here we investigate the ef- 
fect of the window function by estimating the cosmological 
parameters from the VSA dataset with each set of window 
functions. We have used the full compact and extended ar- 
ray dataset. The binning scheme and maximum likelihood 
bandpowers are shown in Table 6. The bandpower window 
functions have been computed using the CBI software (My- 
ers et al. 2003) which has been adapted for use with the VSA 
specifications, including the observing frequency, bandwidth 
and primary beam size. Figure 5 shows the variance and win- 
dow functions for the VSA. Each variance window function 



Table 6. The CMB band-powers (in /iK 2 ) for the complete VSA 
data set combining both compact and extended array data. 



Bin 



Grange 



4ft Tq£(£ + l)C^/27r [fiK 2 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 



100 
190 
250 ■ 
310 
370 
450 
500 - 
580 
640 ■ 
700 
750 
850 - 
950 - 
1050 
1200 
1350 ■ 



190 

- 250 

■ 310 
370 

- 450 

- 500 

- 580 

- 640 
700 

- 750 

■ 850 
950 
1050 

- 1200 
1350 

■ 1700 



156 
220 
281 
333 
410 
475 
537 
611 
670 
721 
794 
902 
987 
1123 
1267 
1440 



3626 
5561 
5131 
2531 
1570 
1811 



1616 
1150 
1561 
1232 
1123 
959 
+438 
411 
+246 
219 
3X3 
356 
2212^6 
1736±356 

1614l|» 
1628l«J 
2486^01 

,+274 
-274 
-+274 
-246 
7+274 
-246 

QOy + 356 

ao ' -329 

75S+ 657 
' 00 -603 



1553Z 
1135 j 
677 4 



Table 7. The priors assumed for the basic parameters. The no- 
tation (a, b) for parameter x denotes a top-hat prior in the range 
a < x < b. 



Basic parameter 



Prior 



e 



ln(10 10 A s ) 



(0.005,0.100) 
(0.01,0.99) 
(0.5,10.0) 
(0.01,0.50) 
(0.5, 1.5) 
(2.7,4.6) 



is normalised to unit area. The bandpower window functions 
are, by definition, normalised to unit area within the band 
limits and to zero outside of the band. The bandpower win- 
dows show negative 'sidelobe' features which indicate the 
anti-correlations of Ces in adjacent bins. 



4.1 Parameter Estimation 

We consider the set of cosmological models described by 
the following six free parameters: the physical baryon den- 
sity LOb — fib ft 2 ; the physical dark matter density u>dm = 
fidmft 2 ; the ratio of the sound horizon to the angular di- 
ameter distance 0; the optical depth to the surface of last 
scattering r; the amplitude of scalar modes A s and the spec- 
tral index of scalar modes n s . We assume spatial flatness, 
setting the curvature density fi k = 1 — Q t ot = 0. We also 
set w — — 1 describing the equation of state of dark en- 
ergy (p — wp) and do not consider tensor modes or massive 
neutrinos in this analysis. 

The parameter estimation was performed using the Oc- 
tober 2004 version of the CosmoMC software package (Lewis 
& Bridle 2002). This uses a new parameterisation with 9 
instead of Ho as a basic parameter. As 6 is less correlated 
with other parameters than H , this has the advantage of 
allowing the Markov chains to converge more quickly. 
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Table 8. Parameter estimates and 68% confidence intervals from 
the marginalised distributions. When the marginalised posterior 
for a parameter does not contain a peak, 95% confidence limits 
are given. 





Variance windows 


Bandpowcr windows 




n 029+ 006 


027+ 006 
u - uz( ^0.006 




0.14±g;gJ 


0.13±8;« 




0.42±g : ?9 


o-39±8:?? 


Q A 


u.oozt 09 


0-61±°i? 


e 


1.04±g;0| 


i nc;+0.02 
l.UOiQ 02 


h 


0.69±g : }8 


n 7-|_|_0.()9 

u -' i=t o.io 


n s 


0.90±8;°; 


qq+0.07 


lnl0 10 A s 






T 


(0.04,0.47) 


(0.04, 0.48) 




21±ig 


22±i} 


Age 


13.1±°;g 


13.2±8;| 



In addition to the priors on the basic parameters, listed 
in Table 7, we also impose priors on some of the derived 
parameters. Specifically, we use top-hat priors on the age of 
the Universe lying between 10 and 20 Gyr, and of Ho lying 
between 40kms~ 1 Mpc~ 1 and 100 kms -1 Mpc -1 . 

The CosmoMC software was run on a 24-node linux clus- 
ter. The chains were run until the largest eigenvalue returned 
by the Gelman-Rubin convergence test reached 0.08. After 
burn-in a total of more than 100,000 samples were collected. 
As successive samples in a Markov chain are, by nature, cor- 
related, the samples were thinned by a factor of 25 resulting 
in approximately 4000 independent samples. These samples 
were then used to calculate the marginalised distributions 
and parameter estimates (see table 8). 

Figure 6 shows the marginalised distributions for each 
parameter. It is clear that no significant bias has been intro- 
duced as a result of using variance windows to approximate 
bandpower windows. The largest discrepancy is seen in the 
parameter estimate for ui b where the bandpower windows 
reduce the estimate by one-third of the 1-sigma error. For 
all other parameters the estimates are consistent to a much 
smaller fraction of the error. The width of each distribution 
is also shown to be unchanged by any significant amount. 
Furthermore, the correlations of the parameter estimates are 
also consistent using both methods. 

These results are perhaps to be expected given the 
sparse nature of the covariance matrix. For the VSA ex- 
tended array data, approximately 5% of the elements of 
the matrix are non-zero and in the limit of a diagonal co- 
variance matrix the bandpower and variance windows are 
equivalent. However, the validity of the approximation also 
depends on the signal-to-noise level and on the binning used. 
The level of correlations between the bins and the gradient 
of the power spectrum across the bin also has an impact 
on the bandpower predicted. For example, using a standard 
ACDM power spectrum and the binning used above, we can 
compare the bandpowers predicted from each window func- 
tion. In the majority of bins, the agreement is within 5% 
but there is a 13% difference in bin 3 where the gradient 
of the power spectrum is steep and there is a strong anti- 
correlation with bin 4. The effect of doubling the size of the 
bins increases the discrepancy between the bandpower pre- 



dictions to a maximum of 28%. Differences of this size would 
be likely to bias the parameter estimates. However, the level 
of agreement between the bandpower predictions is not a 
simple function of bin size and will fluctuate depending on 
the exact binning used. 

The VSA is currently undergoing an upgrade to increase 
the sensitivity of the instrument and to enable a measure- 
ment of the power spectrum up to £ — 2500. The upgrade 
will involve the fitting of larger mirrors with a diameter of 
0.55 m. In this case, the fraction of non-zero elements in the 
covariance matrix will increase slightly to 7.5%. The suit- 
ability of the approximation also depends upon the signal- 
to-noise achieved. The upgraded VSA will have increased 
sensitivity, but taken in conjunction with the reduced signal, 
then the overall signal-to-noise level will remain about the 
same. Given the significant difference between the expected 
bandpower values which may arise and to ensure that this 
does not bias future VSA parameter estimates, analysis of 
the upgraded VSA data will be carried out using bandpower 
window functions. 



5 CONCLUSIONS 

In this paper we have presented three sets of coincident CMB 
observations observed at different epochs by the VSA and 
CBI telescopes. We have chosen to analyse the full datasets, 
rather than focusing on the power spectra, in order to inves- 
tigate any possible systematic effects which may not other- 
wise be revealed. The correlation of the datasets from each 
group was found to be as expected for the 02H and 14H 
mosaics. In the third mosaic the data disagreed with the 
Monte Carlo simulations at a level of 2a. However, this is 
consistent with the phase calibration errors expected from 
VSA data during the summer months. It has been estab- 
lished that this does not affect the power spectrum estima- 
tion due to the correlated nature of the phase errors. The 
results of this analysis reaffirm that both groups have cor- 
rectly characterised the noise properties and systematics of 
the telescopes as well as other, potential data contaminants. 

We have investigated the use of variance windows as an 
approximation to bandpower windows, for the VSA, and 
found that for the data obtained so far, this is a valid 
approximation. We note that alternative binning schemes 
may reduce the suitability of this method and plan to use 
bandpower window functions for parameter estimation from 
super-extended VSA data. 
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Figure 6. The one-dimensional marginalised probability distributions for cosmological parameters estimated using variance windows 
(solid lines) and bandpower windows (dotted lines). 
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